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Abstract 

Using the gauge-gravity duality we have obtained the potential between a heavy quark 
and an antiquark pair, which is moving perpendicular to the direction of orientation, in a 
strongly-coupled supersymmetric hot plasma. For the purpose we work on a metric in the 
gravity side, viz. OKS-BH geometry, whose dual in the gauge theory side runs with the 
energy and hence proves to be a better background for thermal QCD. The potential obtained 
has confining term both in vacuum and in medium, in addition to the Coulomb term alone, 
usually reported in the literature. As the velocity of the pair is increased the screening of the 
potential gets weakened, which may be understood by the decrease of effective temperature 
with the increase of velocity. The crucial observation of our work is that beyond a critical 
separation of the heavy quark pair, the potential develops an imaginary part which is 
nowadays understood to be the main source of dissociation. The imaginary part is found 
to vanish at small r, thus agrees with the perturbative result. Finally we have estimated 
the thermal width for the ground and first excited states and found that non-zero rapidities 
lead to an increase of thermal width. This implies that the moving quarkonia dissociate 
easier than the static ones, which agrees with other calculations. However, the width in our 
case is larger than other calculations due to the presence of confining terms. 


1 Introduction 

The in-medium behaviour of heavy quark bound states is used to probe the state of matter in 
quantum chromodynamics (QCD), where the screening of the potential results in the suppres¬ 
sion of the yields of heavy quarkonium states in relativistic heavy-ion collisions [1]. The heavy 
quarkonium bound states treated nonrelativistically possess the well separated energy scales, viz. 
the heavy quark mass (mg) the momentum transfer (1/r oc mQag) ^ the binding energy 
[E oc rngof) [2]. In addition, the (thermal) medium introduces scales to the previous ones, viz. 
T, the Debye mass, mu (gT) etc. which are also separated (T ::§> mr)) in weak coupling regime 
{g < 1)[3]. The above hierarchies facilitate to develop the sequence of effective theories namely 
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NRQCD, pNRQCD etc. from the underlying theory of strong interaction (QCD) after integrat¬ 
ing out the successive energy scales for both T=0 and T 7 ^ 0[4]. Now if the QQ bound state 
moves with respect to the medium, the synthesis of effective theories becomes more complicated 
due to the additional scales associated with the motion of the pair. However, EFT have recently 
explored the in-medium modifications of heavy quarkonium states moving through a medium for 
two plausible situations: uiq :$>l/r':^T:$>E^mD and mg ':^T:$>l/r,m£)^E [5], which 
are relevant for moderate temperatures and for studying dissociation, respectively. They found 
that the width decreases with the velocity for the former situation whereas for the latter regime 
the width increases monotonically with the velocity. 

The aforesaid hierarchies assume the weak coupling limit and its identification with the rela¬ 
tive velocity. Therefore one needs the framework in strong coupling regime, where the potential 
can be extracted from the Euclidean correlators calculated in the lattice but the limited sets 
of data in addition to the technical difficulties of lattice calculations, limits the reliability of 
the results. Thus some complementary methods for strong coupling are desirable, where the 
AdS/CFT conjecture provides such alternative [ 6 , 7, 8 ] to calculate the potential for a class of 
non-abelian thermal gauge field theories. The expectation value of a particular time-like Wilson 
loop dehnes the potential between a static quark and antiquark at finite temperature. This is 
where the AdS/CFT correspondence makes the calculation easier by mapping the evaluation of 
a Wilson loop in a hot strongly interacting gauge theory plasma onto the much simpler problem 
of Ending the extremal area of a classical string world sheet in a black hole background. The 
first calculation was done by Maldacena [9] for J\f = 4 SYM for T=0 and was later extended to 
finite temperature in [ 10 , 11 ]. 

The physics of quarkonium dissociation in a medium has undergone refinements over the last 
decade [13, 42, 43, 44, ] There are mainly two processes in QCD, which contribute to the ther¬ 
mal width: the first process is the inelastic parton scattering mediated by the space-like gluons, 
known as Landau damping and the second one is the gluon-dissociation process which corre¬ 
sponds to the decay of a color singlet state into a color octet induced by a thermal gluon [14]. 
Recently two of us derived the imaginary component of the potential using thermal field theory, 
where the inclusion of confining string term makes the (magnitude) imaginary component larger, 
compared to the medium-modification of the perturbative term alone [15, 16, 17] However the 
results referred above about the imaginary component of the potential are limited to the weak 
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coupling techniques. Using holographic correspondence, authors in [18, 20] in strong coupling 
limit have obtained the imaginary part of the potential for A/” = 4 supersymmetric gauge theory 
beyond a critical separation, rc of a static QQ pair, by analytically continuing the string configu¬ 
rations into the complex plane. They suggested that the complex-valued saddle points beyond Vc 
{=0.87zh) may be interpreted as the quasi-classical configurations in the classically forbidden re¬ 
gion of string coordinates, analogous to the method of complex trajectories used in quasi-classical 
approximations to quantum mechanics [21]. Complex valued string coordinates at the extrema 
have also been considered in the calculation of total cross section for deep inelastic scattering 
(DIS) in [19]. In an alternative approach in AdS/CFT framework, the imaginary part is calcu¬ 
lated in a strongly coupled plasma through the thermal worldsheet fluctuation method [22]. In 
the aforesaid calculations, the dual geometry taken was the pure AdS black hole metric, hence 
their dual gauge theory is conformally invariant unlike QCD which depends on scale. 

A resurrected interest in the properties of bound states moving in a thermal medium arose in 
recent years due to the advent of high energy heavy-ion colliders. Using holographic correspon¬ 
dence, the potential between two heavy quarks moving through the medium or, equivalently, 
two heavy quarks at rest in a moving medium is calculated in a pure AdS black hole back¬ 
ground [23, 27, 28, 25, 26], where the solution above a critical separation of the pair is abandoned 
and is discussed in details in [24]. Recently the velocity dependence of the imaginary potential 
of the traveling bound states is also calculated [29], following the idea of thermal worldsheet 
fluctuation method [22, 30]. However in abovementioned calculations the metric used is AdS 
black hole, thus the dual gauge theory is conformal, unlike QCD. 

Our aim is to improve the abovementioned calculations in two folds: first of all we improve 
the dual gravity which is somewhat closer to QCD than pure AdS geometry in the sense that it 
accounts for RG flow [31]. Secondly, unlike the world-sheet thermal fluctuations, we delve into 
the solution beyond the critical separation of the pair, which leads to the complex-valued string 
conhguartions and hence the potential turns out to be complex. Recently we have implemented 
both points in a calculation for a static pair [32], and in the present work, we wish to extend the 
calculation for a moving pair. Our paper is organised as follows : In Section 2, we have revisited 
the OKS-BH geometry in brief and in Section 3, we employ this geometry to obtain both the 
real and imaginary parts of the potential by the Nambu-Goto action. Thereafter we calculate 
the thermal width for the ground and excited states in Section 4. Finally we conclude in Section 
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5. 


2 Construction of dual geometry 

In this section we discuss briefly about the OKS-BH geometry that we will use to calculate the 
potential. The discussion presented here closely follows [33, 31, 34] and the readers already aware 
of this geometry may go directly to next section to see the calculations. 

A gauge theory is said to be conformal if it does not flow with the scale, i.e. has a trivial 
renormalisation group (RG) flow, and a conformal theory in a (3 + l)-dimension can be con¬ 
veniently described on the boundary of pure anti-de Sitter space [6]. But if the theory has a 
non-trivial RG flow like QGD, which is confining in IR while conformal in UV, we cannot de¬ 
scribe the full theory on the boundary of some five dimensional space and hence need to look 
differently at running energy scales. One way out is to embed the D branes in the geometry and 
as a result the corresponding gauge theory exhibit logarithmic RG flow. Such a construction was 
done in the Klebanov-Strassler (KS) warped conifold geometry [35]. The gravity dual for KS 
model is warped deformed conifold with three-form type IIB fluxes and the corresponding gauge 
theory is confining in the far IR but is not free at UV. Other problems with KS model include 
the absence of quarks in the fundamental representation and the non applicability of model at 
non-zero temperature. We look at these one by one. 

In stringy model of gauge dynamics, inclusion of fundamental matter is possible by embedding 
a set of flavor branes in addition to the color branes. The strings connecting to the color branes 
in the adjoint representation of U{Nc) group gives gauge particles while those connecting to the 
flavor branes in the adjoint representation of U(iV/) gives the mesons, and those connected to 
both the flavor and color branes in the fundamental representation gives the quarks and anti¬ 
quarks, having the strings with the opposite orientation. 

Theoretically one could go for both large W and Nf in near horizon limit and translates 
the branes into fluxes and derive the gravity background that is a holographic dual of a gauge 
system with gluons and quarks. In practice one does it using probe approximation, where the 
backreaction of the probes on the background can be neglected for the case Nf Nc and the 
flavor physics is then extracted by analyzing the effective action, namely DBI action plus GS 
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action that describes the flavor branes in the color background. This procedure was already 
adopted in AdS^ x model [46], in conflning model [47] etc. 

To introduce the fundamental quarks in the original Klebanov-Strassler model by embedding 
D7 branes in the gravity dual [35, 49, 50] is subtle because the full global solution that incorporate 
back reaction of the D7 branes on the KS background is complicated and so far we only know 
how to insert coincident D7 branes in the Klebanov-Tseytlin background [35]. The resulting 
background is the Ouyang background discussed in [35] which has all the type IIB fluxes switched 
on including the axio-dilaton. In this OKS (Ouyang-Klebanov-Strassler) model, the local metric 
was computed by incorporating the deformations of the seven branes when these branes are 
moved far away from the regime of interest. Since the seven branes are kept far away, the axion- 
dilaton vanish for the background locally. However there will be non-zero axion-dilaton globally. 
The local back reactions on the metric would modify the warp factors to the full global scenario. 


Now to switch on a non-zero temperature we need to insert a black hole into this background 
and the Hawking temperature will correspond to the gauge theory temperature. Combining all 
physics ingredients together the metric in OKS-BH geometry looks like 




Vh 


+ Vh 




dMl 


( 1 ) 


gi{u)dt^ -|- dx^ + dy^ + dz^ 

where h is the warp factor, gi{u) are the black-hole factors, u denotes the extra dimension and 
dA4| is the metric of warped resolved-deformed conifold. 


The above picture works well at IR. In fact OKS model has emerged as the most promising 
candidate to study large N strongly coupled QCD, but the effective number of degrees of freedom 
grow indefinitely in the UV. Moreover in the presence of fundamental flavors such a proliferation 
leads to Landau poles and hence the UV divergences of the Wilson loops. Our aim here is to 
detemine the Nambu-Goto action with a dual gravity of a gauge theory at finite temperature, 
which resembles the main features of strongly coupled QCD, i.e., a gauge theory which is almost 
conformal in the UV with no Landau poles or UV divergences of the Wilson loops, but has 
logarithmic running of coupling in the IR. In other words we want AdS-Schwarzschild geometry 
in asymptotic limit and therefore we need to add the appropriate UV cap to it. 


Now the addition of UV caps in general can change IR geometries. However the far IR 
geometries have not been changed by the addition of UV caps because the UV caps correspond 
to adding non-trivial irrelevant operators in the dual gauge theory. These operators keep far IR 
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physics completely unchanged, but physics at not-so-small energies may change a bit. So the 
IR geometry part has been modihed to obtain the desired dual gauge theory in the following 
steps [31, 33, 34]: 

i) switch off the three-form fluxes at large r so as to keep only five-form fluxes to get the AdS 
type solutions. 

ii) modulate the axio-dilaton behaviour so that the axio-dilaton held decay rapidly at large r 
instead of the logarithmic growth, to get rid off the Landau poles. 

hi) implement both of them by switching on anti D5-branes and anti seven-branes with electric 
and magnetic huxes on them to kill off the unnecessary tachyons and thereby restore the stability. 
Additionally the anti seven-branes, along with the seven-branes, were embedded in a non-trivial 
way so as not to spoil the small r, or the IR, behaviour of the theory. The anti D5-branes 
are dissolved on the seven-branes, but they contribute to the three-form fluxes. Both of them 
together control both the three form fluxes as well as the axio-dilaton at large r. 


As a consequence, the above metric will receive corrections as 


ds^ = 


Vh - 


g{u)dt^ -I- dx'^ -I- dy'^ -|- dz"^ 


Vh 


17 (^) Quudu “t“ Qrnndx dx 


( 2 ) 


where we have set black hole factors gi{u) = g2{u) = g{u). The corrections guu are of the form 
and appear in the metric because the existence of axio-dilaton and the seven-brane sources 
tell us that the unwarped metric may not remain Ricci flat. Thus the corrections may be written 
as 

OO 


9uu - 




^uu^i 


( 3 ) 


i=0 

where a^u,! are the coefficients independent of the coordinate u and can be solved for exactly as 
shown in [31]. The warp factor, h can be obtained as 


h = 



OO 


CLi 


i=l 

where Oj are coefficients, again independent of the coordinate u. These are of 0{gsNf) and can 
be solved for exactly as shown in [31] and L denotes the curvature of space. This metric reduces 
to OKS-BH in IR and is asymptotically AdS^ x M 5 in the UV. It describes the geometry all the 
way from the IR to the UV. 
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With the change of coordinates z = 1/u, we can rewrite the metric (2) as 

ds^ = = Anz'^-^ [-g{z)df + d'^^] + - —(4) 

where is the metric of the internal space and A„’s are the coefficients that can be extracted 
from the Oj’s as follows: 


1 

7 ^ 


L'^z’^y/ttiZ^ 


— An.Z 


n—2 


1 


aiz 

ao-^ + 




( 5 ) 


which gives Aq = j%,Ai = —■^,A 2 = so on. Note that since a^s for i > 1 

are of 0{gsNf) and oc y/ggN, so in the limit gsNj —?■ 0 and —?■ cxd all Aj’s for i > 1 are very 
small. The corrections mentioned above in Eq.(2) are accommodated via Biz^ series which 
is given by 

Biz^ = 1 + ( 6 ) 


In fact the complete pictnre can be divided into three regions. Region 1 is the IR region 
where we have pnre OKS-BH geometry. Region 3 is the UV region where UV cap has been 
added. And region 2 is the interpolating region between UV and IR. The background for all 
these three regions and the process of adding UV cap has been described in full details in [31]. 
Also the RG flow associated with these regions and the corresponding field theory realizations 
have been discussed in [36]. We shall not go into the complete details here and will use the metric 
given in Eq.(4) in extremizing the action to calculate the potential in the next section. 

Another model that is extensively used to study certain IR dynamics of QCD is the Sakai- 
Sugimoto model[48] which is a beautiful brane construction in type IIA theory. The model 
consists of a set of N wrapped color D4-branes on the circle, while the flavor branes D8 and 
D8 placed at the anti-nodal points of the circle give the mesonic states. In the gravity dual, 
the wrapped D4-branes are replaced by a geometry, i.e an asymptotically AdS space, but the 
eight-branes remain and so does the circular direction. However the Sakai-Sugimoto model does 
not have a UV completion, which is very important for the study of heavy quarkonium. The UV 
contribution to the the heavy quark potential for the heavy quark bound states are dominant 
compared to IR part because since their mass is larger than QCD scale! so the relative velocity 
is very small and becomes weakly coupled {v ~ as, as is the strong interaction coupling). The 
comparison of Sakai-Sugimoto model with the one that we use here has recently been done in 
151], 
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3 Potential between a heavy quark and a heavy antiquark 


The potential at finite temperature can be obtained from the Wilson loop’s expectation value 
evaluated in a thermal state of the gauge theory : 

lim {W{C)) ~ (7) 

T—>oo 

According to the gauge/gravity prescription [ 6 ], the expectation value of W{C) in a strongly 
coupled gauge theory dual to a theory of gravity is 

( 8 ) 


where Sgtr is the classical string action propagating in the bulk evaluated at an extremum. Hence 
the potential can be obtained by extremizing the world-sheet of open string attached to the heavy 
quark pair located at the boundary of AdSs space in the background of AdSs space and in the 
background of AdSs black-hole metric, for determining the potentials at T = 0 and T 7 ^ 0, 
respectively. 

We consider a quark-antiquark pair [QQ) moving along direction with some velocity v. It 
is convenient to boost the system to the rest frame (T, Xg) of the quark-antiquark pair as 

dt = dt' cosh 7] — dx'^ sinh r], 

dx3 = —dt' sinh 7] -\- dx'^ cosh r], (9) 


with tanh rj = v. The QQ dipole is now at rest but the quark- gluon plasma is moving with 
velocity v in the negative x'^- direction. Thus under the boost, the metric (4) reduces to 


ds^ = AnZ 


n—2 


dt ( —1 H— 4 cosh r] 1 -|- dx{ -|- dx^ + dx^ ( 1 - 1 — j sinh rj 


— 2 dt dxs sinh 77 coshr; 


+ 


Biz^ 


AmZ"^+‘^g{z 


-dz^ + 




( 10 ) 


Now the Nambu-Goto action can be defined in terms of string coordinates {a, r): 


Sng = ^ y dadr^-det , (11) 

where 0 is the background dilaton field, which is responsible for breaking of conformal symmetry 
of theory and is given by 


<!> = logs, - 9,£>„+m.z”+“" 


( 12 ) 







We parametrize the two-dimensional world sheet as 

T r d d 

T = t, a = e X 2 = const, X 3 = const, z = z{xi), da = Try db = — , (13) 

22 OT da 

thus the Nambu-Goto action can be rewritten as 


r 


Sng = ttz I dxi\l{AnZ^y ( — - ^cosh rj ) + BmZ' 


1 1 


271 ' \z^ zt 7 ( 7-^7 


{z'y (7 — z"^ cosh^ rj) 


where we have used 


= Biz' + 2g{z)gl(n + . 


(14) 


The Nambu-Goto action can also be written as an integral over z. 


T dz I / 1 1 i 2 \ n. (^ 0 ^ (7 - 2 ;^ cosh^ 77 ) 

Sng = —I —\l(AnZ^)^{— -7cosh^?71 -h 1 -r-i -- (15) 


TT 




(7 - ^7 


Since this action does not depend explicitly on x, so the corresponding Hamiltonian will be a 
constant of motion, i.e., 

fjT 

H = z'—-L = Co (say) , (16) 

where the constant Cq can be determined from the following equation: 


1 f,j^ 2 BmZ^ ( 7 -7 cosher/) ^2 


(7 - ^7 


1 1 


1 = - ( ^ ^ cosh^ 7] 1 ^^"7 ^ = Co 


(AnZ 


n\2 


Since at z = Zmax, z'=0, where z^ax is the maximum of the string coordinate along the hfth 
dimension, so we can hnd out the derivative z' by the simple algebra 


z = 


dz 

dx 


{AnZ^f Z, 


n\2 ^2 

max 


(7-7)(7-7cosher/) 


\/BmZ'^ Zl y {zl - Zf^^^ cosh 7]) 


X 


i^nZ^ax) {4 - ^rnax COsh r/)7 


. 4 \ 2 


''max 


{AnZ'^y{zl^ — 7 cosh^ r])zfi 
Integrating both sides, the above equation yields the separation of the QQ pair, r as 


r = 


X 


27\/7-4ax cosher] (H„7 


(4 - ^tax cosh^ h)(^nCaJ7^ 

(7 - 7 cosh^ r]) {AnZ^y z^^^ 


dz 


-4K 


(17) 


\/( 7 -^ 7 ( 7 -^^ cosher/) (AnZ^y 


(18) 
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Expanding the sqnare root and keeping only the first term because the coefficients Aj’s are small, 
the separation (r) becomes 


r = 


2zl{AnZl^ax) \/ COsll^ T] 




where the integral I is defined by 


I = dz 


■VK, 


(4 - “ ^‘‘cosh^r/) 


(19) 


( 20 ) 


Now we obtain the action 


Sng — 



dz IBmZ"^{zf^ — z^ cosh^ r]) 

(4^) 


A _ {4-^Lx cosher/) 

V {4 - cosh^ rj) {Ar^z^yz^^^ 


1 

2 


Expanding the square root and keeping only the first two terms (since the coefficients Aj’s are 
small), the above action can be written as 


5 

o/vn — — 


TT 



^""(4 - cosh^ r/) l(zl- 4ax cosh^ 




, ( 21 ) 


where I is the same integral as in (20). Thus substituting the integral I in terms of r, the action 
becomes 




NG 


T_ 

TT 


dz BmZ'^{z'l — z^ cosh^ rf) 


T \l4-^max cosher/ 
47r 


2^2 

'^h max 


4 r 

^n^max ’ 


= S^ + S^ 


( 22 ) 


The hrst term in the action, diverges in the lower limit of the integration, so we regularize it 
by integrating from e (instead of 0) to z^ax [31] and then identifying the divergent term in the 
integral 

Zmax 


Let us also assume for simplicity that the higher coefficients Ai and Bi are very small for i > 3 
and can be neglected. This simplification reduces the series AnZ"' to 1 + A 2 Z^ and the other series 
BmZ'^ to 1 + B 2 Z^, where we take Aq = 1 and Ai = 0 and similarly Bq = 1 and Bi = 0. This 



dz\fW„ 


' zf^ — z'^ cosh^ 
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will also simplify all the expressions and helps us to keep an analytic control on the equations. 
Therefore, in the limit of small coefhcients, becomes 


S^ = - 

71 


dz{l + zt — z^ cosh^ T] _ T 

= 7 


4 - 


dz F{z, v) , 


where the integral, F{z, v) can be expanded in the Taylor series in velocity: 

r./ ^ r./ X dF Id'^F 2 

F(z,v) ^ F(z,0) + V + . 


(23) 


We can omit the dotted terms in the small velocity limit and can write as 

,2 pZmax / D \ 1 

-02 ^ \ 1 

T 


= - 


TT 




dz{z‘^ + 




(24) 


After subtracting the divergent piece in the limit e —)■ 0, we get the renormalised Nambu-Goto 
action 

1 


T 

nren _ ' 

TT 


B2 

^max 

2 


^(-4B2Z,z^a. + 2(-2 + i?2^^)tan-' 

IGZh V \ Zh 


(2 + B 2 Zf,) log 


2\1 {^h + ^max ) 


{zh — zmax) 
and hence the potential is obtained by 

cren 


+ 


T \l^h-^max cosher; 

47r 


2p‘ 7^ 

max 


A 7^ r 
^n^max ’ 


(25) 


Vqq = lim 
«« r^oo T 


Bo 


^ '^max 
2 Zn 


/ 


4:B2ZhZmax + 2(-2 + 52^) tan' 


— 1 / ^max 


(2 + B 2 Zf^) log 


2 m {^h + ^max) 


{zh — zmax) 


+ 


1 \l^h-^^ax cosher/ 


Zh 


47r 


max 


AnzT^ r . 


(26) 


The potential thus obtained are functions of both Zy^ax and r, which are redundant. To make 
the potential as a function of r only, we will express Z max as a function of r and then plug in to 
the above potential. To do that we hrst concentrate on the integral I (20) which, keeping upto 
the second-order in both series in z, reduces to: 

_^Vl -b B2Z'^ _ 

{z^ — z'^){zf^ — z‘^ cosh^ r]) {1 + A 2 z‘^)'^ 

^max 

dz r{z,v) (27) 
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In the limit of small velocity, we expand the integrand I'{z,v) in the Taylor series in velocity, v 


as: 

w/ ^ ru ^ dr id^r 2 
and neglecting the higher order terms, the integral can be written as: 

, z^\/\ + r^max + B2Z^ 

Wr_ -I- _ I (12^ _ 

(4 - ^4)(1 + A^zY 2 io (4 - + ^2^2)2 

We will now hnd the solution for the separation, r in two limits, namely Zmax » Zh and 
Zmax « Zh- First for the Zmax » Zh limit, the separation r from Eq.(19) becomes as a function 
of z 

^max 





2tzl 


Boil - —) 

= “(ci + iC2 - —(C3 + iC4)){A2z‘^^^ + 1) + ^ 

/ 0/12 Z^ax 


+ 


1 zt A,_{1 - v^), v\ , . 

2 -“ ¥<‘=3 + +- sn ^2 ,3 - 

^ '^max ^ 

(l-n^) zi 


^2 '^max 


(ci + iC 2 - 17 (C 3 + iC4)) 


where Ci , C 2 , C 3 and C 4 are dehned as 


71 

Cl = w 


1 { 2 -B 2 ZI) {QA 2 -B 2 + Alzt{2A2-3B2)) 


C2 = 


C 3 = 


C 4 = 


8 [zh {A 2 ZI - 1)2 
TT {2 + B2ZI) 
8Zh {A 2 ZI + 1)2 


v^(-i + ziiz^y 


TT 


{ 6 A 2 -B 2 + Alzt{ 10 A 2 - 7 B 2 )) (6 - 5B2zI + 242Z^(2 + B 2 ZI)) 


V^(-l + 

7r(6 + 5-02^^ + 2422:^(—2 + i?2^^)) 


Az,iA2zl - 1)3 


32zh{A2zl + 1)3 
We now invert the series (28) to obtain z^ax in terms of r as. 


\J 2 A 2 z‘f^{—ici + C 2 + \{icz — C 4 )) 

(-ici + C 2 + y(*C3 - C4))zl 


i B 2 zl{2 - v^) 


2 A 2 . 


A^(l-v-^) 12 2l2 ry(l-t,2) ’ 


(28) 

(29) 

(30) 

(31) 

(32) 


(33) 


which shows that the string coordinates become complex^ and hence the potential from Eq. (26) 

^ This solution was earlier abandoned in evaluating the potential for the heavy quark pair in a moving 
medium [24]. 
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Figure 1: Variations of the real and imaginary part of the potential with the separation r at various velocities 
(at a constant temperature 400 MeV), where both coefficients, A 2 and B 2 are taken as 0.24. 


becomes imaginary: 


^QQ = -[^Zmax{l +^) - + B2Z^-i{2 + B2Zf^)) - - -(^ + y) “ TF3— 

n \ Zj Zj ^max ^ ^max 


ir 


^2Zmax + 1 


A 2 zi{l-v^) Iziil-v^) 


(34) 


After substituting Zmax in terms of r from Eq.(33), we obtain the real and imaginary parts of 
potential as 


Re[V{r,T,v)] = 


r‘^ (ci - yCs) 


87r2;^( (C2 - ^^ 4 )^ + (ci - 


+ 


(10 — — 6n^) 

24(1 -n2) 


X 


\ 


(c2 — Y'^ 4 )^ + (ci — (^2 ~ 


Bo 


2 ( (C2 - ^^4)^ + (ci - f-ca)^) 


V^^TlZh 


’\/l — 


16zh 


[-2 + B 2 ZI) - 


(4^2 + B2)zh (18 - - lOn^) 


\/2A^7r 


40(1 -n2) 


X 




(c2 ~ ~ (^2 ~ 


(35) 


rv i — n 


and 


Im[V{r,T,v)] = 


(c2 - Y^y 


Snzf^i (c2 - ^€ 4)2 + (ci - n 2 c 3 ) 2 ) 


+ 


(10 — — 6n^) 

24(1 -n2) 


X 




(C2 - Y^ 4 .y + (ci - Y^yy - (C2 - Y^^) B2 


2 ( (C2 - ^C 4)2 + (Ci - ^C 3 ) 2 ) 


y/lA^TTZh 


ry/l — ^2 (36) 
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Figure 2: Variations of the thermal width for the charmonium ground state with the temperatures at increasing 
velocities (left) and with the velocity at increasing temperatures (right). 


respectively. We have plotted the variation of the real and imaginary parts of potential with the 
separation r for different velocities. We hnd that the potential becomes stronger with the increase 
of velocity, which can be understood from the fact that the effective temperature gets reduced 
due to the recession of the medium and hence the screening becomes weaker. The imaginary part 
increases in magnitude as the velocity increases. We will now hnd the potential for the other 
extreme limit he., very small Zmax and Zmax « Zh- In this limit, the separation r becomes 


r = 


max 


( 3 i ?2 ~ 2^42) 
15 




3 

max 5 


(37) 


which will eventually be inverted to express Zmax in terms of r as 


Zmax 


3r 

y 


27 

-(2A^ - 3B,y + 0|r]= 


(38) 


Thus the potential in this limit reduces to 


&m,yr,T)| —V'* - A + £(T_^, + o[r|’ (39) 

Hence upto hrst order in r, the potential is independent of velocity in this limit. The hrst term 
here is like the coulombic term while the second term is the linear conhning term. 
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Figure 3: Variations of the thermal width for the charmonium first excited (2S) state with the temperatures at 
increasing velocities (left) and with the velocity at increasing temperatures (right). 

4 Calculation of thermal width 


The thermal width Fqq arises due to the imaginary component of the potential and is dehned 
as 


rgg = -(c/’llm VQQ(r,T,ti)|V>), 


(40) 


where the wave function for the ground state (IS) is given by 


and for the hrst excited state 


V' = 


p-r/ao 




3/2 


'll) 




(41) 


(42) 


in a Coulomb-like potential, with oq = (p-e^) ^ being the Bohr radius (p is the reduced mass 
and is the square of the electric charge). Even though the real part of the potential is not 
purely Coulombic but the leading contribution for the potential for the deeply bound QQ states 
in a conformal plasma and thus justihes the use of Coulomb-like wave functions to determine the 
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Figure 4: Same as Figure 2 but for the bottomonium state. 


width. Thus we calculate the thermal width for the ground-state(lS) as 


rQQ(i^) = 


drr'^e Im l/QQ(r, T) 


^0 -'0 


9727r'^T^ 


(C2 - f C4) 


200m^ ((C 2 - f C4)2 + (ci - 4 c3)2) 32v'5^^ 


45\/27rTi?2 . 2 1 (10 — — 6u^) 

- (1 - - - 


24(1 -u2) 


X 




+ (ci - ~ (^2 - 

2( (C2 - 4c 4)2 + (Ci - ^C3)2) 


and for the hrst excited state (2S), 
eSOdTT^T^ (C2 


rQQ(2^) = 


C4) 


lOOm^ ((C2 


1035V27rT52 
-C4)2 + (ci-f^ “ 2567Spr 


(1 


,2'i4 


(43) 


i (10 — — Ou"^) 


24(1 


X 


\ 


+ (ci - ^C3)2) - (C 2 - YC 4 ) 
2( (c2 — ^^4)^ + (ci ~ 


(44) 

2 I \'^i 2 j 

where the second term is the higher-order correction in velocity and contribute smaller to the 
hrst term in the limit of small velocities (which is also checked numerically) and thus can be 
neglected. In the limit of small velocity, the thermal width for IS state can be simplihed into 


rQa(i^) 


C2 


9727r5T^ 


200m^ [{4 + 4) 


+ 


C 4 2 c 2 (c 2 C 4 -1- CiC 3 )\ 


(cf + ci) 


(ci + ci)2 ; 


(45) 


We have plotted the thermal width for the ground and hrst excited states of charmonium (in 
Fig. (2) and (3) ), where we found that the width is always broadened with the the increase of 
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Figure 5: Same as Figure 3 but for the bottomonium state. 

temperature and also increases slowly with the velocity. It can be understood from the fact 
that as the velocity of the medium increases, the flux of plasma passing through the QQ pair 
increases and hence causes to broaden the thermal width. The same reason applies to the increase 
of imaginary part of potential with the velocity. As a consequence, the pair will be dissolved 
earlier into a moving medium compared to the static medium. Another observation is that the 
thermal width for bottomonium states (Fig. (4) and (5)) is much smaller than the charmonium 
states because bottomonium states are bound tighter than charmonium states. 

5 Conclusions 

In summary, we have obtained the inter-quark potential at hnite temperature in a dual gravity 
closer to thermal QCD for a heavy quark and antiquark pair moving normal to its orientation. 
When the (hanging) string lying on the hfth dimension, z^ax is far above the horizon, ie. Zmax << 
Zh, the Nambu-Goto action gives rise to the similar form of Cornell potential, unlike the Coulomb 
term alone usually reported in the literature. On the other hand when the string reaches deep 
into the horizon (z^ax >> ^h), the potential develops an imaginary component. Alternatively 
it can be stated that beyond a critical separation of QQ pair, the string coordinates become 
imaginary and as a result, the potential becomes complex. The imaginary part vanishes in the 
limit of small separation which is also seen in the perturbative calculations. Furthermore as the 
pair starts moving, the screening of the potential becomes smaller, which may be understood 
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qualitatively by the decrease of effective temperature. However for a particular velocity of the 
pair, the screening becomes stronger with the increase of the temperature as observed in the 
potential for a static pair. 

Since the quarkonium dissociation is presently thought mainly due to Landau damping in¬ 
duced thermal width, we have therefore obtained the thermal width for the ground and hrst 
excited states from the imaginary part of the potential. We found that the width not only in¬ 
creases with the increase of temperature, it also increases slowly with the velocity of the pair, 
which makes sense. The effect of enhancement of width is more pronounced to the charmonium 
states than the bottomonium states. Since the study of quarkonium bound states propagating 
in the QGP at hnite velocity is nowadays relevant in heavy-ion collisions so our study gives an 
theoretical input in this regard. Indeed the PHENIX Collaboration [37] had found a substantial 
elliptic flow for heavy-flavor electrons, indicating a signihcant damping of heavy quarks while 
travelling across the medium [38]. ^ 
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